Prl-1 

© 

vo" 
o 
o 

£ 

(N ! Non-equilibrium relaxation of an elastic string in random media 
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Abstract 

We study the relaxation of an elastic string in a two dimensional pinning landscape using 
Langevin dynamics simulations. The relaxation of a line, initially flat, is characterized by 
a growing length, L(t), separating the equilibrated short length scales from the flat long 
distance geometry that keep memory of the initial condition. We find that, in the long 
time limit, L(t) has a non-algebraic growth, consistent with thermally activated jumps 
over barriers with power law scaling, U (L) ~ L e . 



Understanding the dynamics of elastic manifolds in random media has been the focus of 
intense activities both on the theoretical and experimental side [1]. Because of the competition 
between disorder and elasticity, glassy properties arise leading to divergent barriers separating 
metastable states. The creep, i.e. the steady-state response to a small external force, has been 
studied theoretically [2,3] and observed experimentally [4]. Much less is known about the non- 
stationary relaxation to equilibrium. Theoretical attempts to tackle this problem have been 
made using mean field and renormalization group approaches [5,6], while numerical simulations 
C$ 1 have been mostly restricted to 2— dimensional systems [7] . Direct applications of these results 
£h . to one dimensional domain walls are however difficult. 

In this paper we study, numerically the slow non-stationary motion of an elastic string 
' ^ , in a two dimensional random medium, relaxing from a flat initial configuration. The string 
q obeys the equation of motion jdtu^z, t) — cd^u(z, t) + F p (u, z) + n(z, t), where 7 is the friction 
£j , coefficient and c the elastic constant, F p (u, z) the pinning force for a random bond disorder 
^ ■ and T)(z, i) is the thermal noise. The details of the numerical method are given in [8]. 

Fig.^a) shows the typical relaxation of a string and Fig.^b) the evolution of the structure 
factor S(q,t) which characterize the geometry of the line [3,8]. At short length scales the 
line has reached equilibrium and it is characterized by the well known roughness exponent 
C = 2/3 [9]. At large length scales a plateau still keeps memory of the initial flat condition. 
A unique crossover length L(t) can be defined. Its evolution is shown in Fig.^c). Two main 
scenarios have been proposed to describe the long time growth of L[t). The first one relies 
on phcnomcnological scaling arguments, based on creep. At low temperatures the relaxation 
is dominated by the energy barriers U(L) that must be overcomed in order to equilibrate 
the system up to a length scale L. Using the Arrhenius thermal activation law we can thus 
express the relaxation time t(L) ~ exp[/3E/(L)]. Even if the exact numerical determination 
of U(L) is an NP-complete problem it is usually conjectured that the typical barriers of the 
energy landscape scale, asymptotically with L, the same way as the free energy fluctuations: 
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Fig. 1. (a) Typical configurations of the relaxing string (b) S(q,t) at T — 0.5. The dashed line 
corresponds thermal equilibrium solution, (c) Characteristic length scale L(t) of a string of size 
L — 256. (o) symbols and the dotted line correspond to the numerical and analytical L(t) for the 
clean system, and (•) symbols to the disordered case. The continuum line fits to the logarithmic 
growth, and the dashed line is a fit to a power law growth at intermediate scales. Inset: exponent 6 
of the logarithmic law extracted from the fit in the time interval ti < t < 10 6 , for L — 512 (A) and 
L = 256 (o). 



U(L) ~ L e , with 9 = 1/3 for a line. Numerical calculations [11] and FRG calculations [2] seem 
to confirm this conjecture. Following these arguments we infer that L(t) ~ [log(t)]» [8,10]. 
The second scenario predicts L(t) ~ i 1 / 2 , like for the clean system but with a new temperature 
dependent exponent z > 2 [6,7]. This power law scaling can be ruled out from our data, due 
the observed bending in the log-log scale in Fig.^c). On the other hand a two-parameters fit 
to the logarithmic growth law gives an exponent 9 which, at long times, becomes size and time 
independent, as shown in the inset of Fig. ^ However we found 9 w 0.49. This value, bigger 
than the expected 9 = 1/3, indicates either a violation of the expected scaling of barriers or 
the presence of logarithmic corrections in the leading term: U(L) ~ L 1 / 3 log M (L) [10, 11]. 
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